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Abstract 

This paper deals with stationary Gibbsian point processes on the plane with 
an interaction that depends on the tiles of the Delaunay triangulation of points via 
a bounded triangle potential. It is shown that the class of these Gibbs processes 
includes all minimisers of the associated free energy density and is therefore 
nonempty. Gonversely, each such Gibbs process minimises the free energy density, 
provided the potential satisfies a weak long-range assumption. 
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1 Introduction 

It is well-known that stationary renewal processes with a reasonable spacing distribution 
can be characterised as Gibbs processes for an interaction between nearest-neighbour 
pairs of points Section 6]. Here we consider an analogue in two dimensions, viz. 
Gibbsian point processes on with an interaction depending on nearest-neighbour 
triples of points, where the nearest-neighbour triples are defined in terms of the Delaunay 
triangulation. Recall that the Delaunay triangulation is dual to the Voronoi tessellation, 
in the sense that two points are connected by a Delaunay edge if and only if their Voronoi 
cells have a common edge. Since the Voronoi cell of a point consists of the part of space 
that is closer to this point than to any other point, this means that the Delaunay 
graph defines a natural nearest-neighbour structure between the points. (Of course, 
the analogy with renewal processes does not reach too far because the independence of 
spacings under the Palm distribution, which is characteristic of one-dimensional renewal 
processes, is lost in two dimensions due to the geometric constraints.) 

There is a principal difference between the Delaunay interactions considered here and 
the pair interactions that are common in Statistical Physics. Namely, suppose a point 
configuration uj is augmented by a new particle at x. In the case of pair interactions, 
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X is subject to some additional interaction with the particles in iv, but the interaction 
between the particles of u is not affected by x. In the Delaunay case, however, the 
particle at x not only gives rise to some new tiles of the Delaunay triangulation, but 
also destroys some other tiles that were present in the triangulation of uj. This so-called 
non-hereditary nature of the Delaunay triangulation blurs the usual distinction between 
attractive and repulsive interactions and makes it difficult to use a local characterisation 
of Gibbs measures in terms of their Campbell measures and Papangelou intensities. Such 
a local approach to the existence of Gibbs measures for Delaunay interactions was used 
in the previous work [21 [3l El [7] and made it necessary to impose geometric constraints 
on the interaction by removing triangles with small angles or large circumcircles. 

In this paper we address the existence problem from a global point of view, which is 
based on stationarity and thermodynamic quantities such as pressure and free energy 
density. Specifically, we show that all minimisers of the free energy density are Gibbsian, 
which implies the existence of Delaunay-Gibbs measures because the entropy density 
has compact level sets. The converse part of this variational principle is harder and 
requires the comparison of different boundary conditions in the thermodynamic limit. 
It is here that the non-hereditary nature of the interaction shows up again, but it can 
be controlled with the help of stationarity and an additional condition which is much 
weaker than the geometric constraints mentioned above. In contrast to [7], however, 
we need to assume throughout that the interaction potential is bounded, and therefore 
do not cover hard-core interactions that forbid particular shapes of the tiles. We note, 
however, that some ideas of the present paper can be used to establish the existence 
of Delaunay-Gibbs measures in a more general setting that includes also the hard-core 
case, see [8]. The extension of the variational principle to such interactions is left to 
future work. As a final comment, let us emphasise that we make repeated use of Euler's 
polyhedral formula and the resulting linear complexity of the Delaunay triangulations, 
and are therefore limited to two dimensions, as was already the case in the previous 
papers mentioned above. 

2 Preliminaries 

2.1 Configurations and Delaunay triangulations 

A subset w of M? is called locally finite if card (a; n A) < oo for all bounded A C M^; 
each such uj is called a configuration. We write for the set of all configurations to. 
The configuration space 17 is equipped with the cr-algebra that is generated by the 
counting variables N/^ : to card{uj D A), with A an arbitrary bounded Borel subset 
of M^. 

For each Borel set A C we write Qa = {w € : w C A} for the set of configu- 
rations in A, pr^ : uj uj\ := lj n A for the projection from 17 to I^a, T'^ = J-\^\ for 
the trace cr-algebra of J- on 17 a, and J-\ = pr^^.Fj!^ for the cr-algebra of events in 17 that 
happen in A only. 

For each configuration cj G 17 we consider the Delaunay triangulation D(w) associated 
to UJ. By definition, 

(2.1) D(a;) = {r C w : cardr = 3, w n B{t) = 0} , 
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where B{t) is the unique open disc with r C dB{T). D(a;) is uniquely defined and 
determines a triangulation of the convex hull of u whenever u; is in general circular 
position, in that no four points of u: lie on a circle that contains no further points of 
oj inside [IT]. If this is not the case, one can apply some determistic rule to make the 
Delaunay triangulation unique. Indeed, let 

T := {r C : cardr = 3} = {{x,y,z) € (M^)^ : x ^ y z] 

be the set of all triangles (or tiles) in where stands for the lexicographic order in 
M^. The triangles in T can be compared by the lexicographic order of (M^)'^, and this 
in turn induces a lexicographic order on finite collections of triangles. Now, if n > 4 
points of Lo lie on a circle with no points inside then the associated Delaunay cell is a 
convex polygon having these n points as vertices. To define a unique triangulation of 
this polygon one can then simply take the smallest among all possible triangulations. 
Confiicts with other possible polygons cannot arise because the tessellations inside and 
outside a fixed convex polygon do not depend on each other. 

Let us note that the prescription lv ^i^^) is a mapping from 0, to the set 0,(T) 
of all locally finite subsets of T. If 0,(T) is equipped with the cr-algebra ^{T) that is 
defined in analogy to one can easily check that this mapping is measurable. 

Next we assign to each tile r G T a centre and a radius. Specifically, for every r E T 
we write c(r) for the centre and q{t) for the radius of the circumscribed disc B{t). 
The centres allow us to consider D(u;) as a germ-grain system, i.e., as a marked point 
configuration of germs in and marks in the space 

To = {r G T : c(r) = 0} 

of centred tiles, by considering the mapping 

(2.2) D:uj^ {{c{t),t-c{t)):t ^D{uj)] 

from J7 to the point configurations on x Tq. Here we write r — c(r) := {y — c(r) : y G r} 
for the shifted tile. 

A crucial fact we need in the following is the linear complexity of Delaunay trian- 
gulations, which is expressed in the following lemma. This result follows directly from 
Euler's polyhedral formula, and is the main reason why we need to confine ourselves to 
two spatial dimensions; see Chapter 11, and [TT], Remark 2.1.4. 

Lemma 2.1 For a simple planar graph on n vertices, the number of edges is at most 
3n — 6, and the number of inner faces is at most 2n — 5. In particular, every triangulation 
with n nodes consists of2n — 2 — d triangles, where d is the number of nodes ( or: number 
of edges) along the outer boundary. 

2.2 Stationary point processes and their tile distribution 

Let ^0 be the set of all probability measures P on (i7, T) that satisfy the following 
two properties: 

(S) P is stationary, that is, P is invariant under the shift group Q = {'&x)xeS^ on Q, 
which is defined hy '■ ^ ^ ^ — x := {y — x : y & uj}. 
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(I) P has a finite intensity z{P) = |A|^^ j N^dP < oo. Here, A is any bounded 
Borel set in M? (which can be arbitrarily chosen due to stationarity) , and |A| its 
Lebesgue measure. 

Each P S is called a stationary point process on with finite intensity. For 
A C M^, we write Pa '■= P o pr^^ for the projection image of P on (Q\,J^j^), which can 
of course be identified with the restriction P\J-'\ of P to the events in A. 

Every P S defines a germ-grain model P, namely the distribution of P under the 
mapping D defined in (j2.2p . That is, P is a stationary marked point process on M? with 
mark space Tq. Let P^ be the associated Palm measure on %xn and fj.p = P°( • X Q) 
the associated mark distribution, or centred tile distribution, on Tq. By definition, 

(2.3) fdx f f,p{dT)f{x,T) = I P{du) f{c{T),r-c{r)) 

for all nonnegative measurable functions / on X Tq. For each P E B^q, ^p has total 
mass II^pII = 2z{P), as follows from Euler's polyhedral formula; see, for example, [n\ 
Eq. (3.2.11)] or |20l Theorem 10.6.1(b)]. 

Let us say a measure P S is tempered if 

(2.4) j \B{T)\fip{dT) < oo 

We write for the set of all tempered P G ^q. Of course, ()2.4p is equivalent to the 
condition J q{t)'^ jipidr) < oo. Moreover, (|2.3p implies that 

(2.5) J \B{r)\f,p{dT) = J dxj Mdr) l{.^Bir)} 

= [ Yl ^{c(r)eB(r-cir))} P{du^) = j cardjr G : G P(t)} P{du) . 

So, P is tempered if and only if the last expression is finite. A sufficient condition for 
temperedness will be given in Proposition 14.91 

The most prominent members of are the Poisson point processes, which will 
take the role of reference processes for the models we consider. Recall that the Poisson 
point process with intensity z > is characterised by the following two properties: 

(PI) For every bounded Borel set A, the counting variable is Poisson distributed 
with parameter 2:|A|. 

(P2) Conditional on A'^a = the n points in A are independent with uniform distri- 
bution on A, for every bounded Borel set A and each integer n. 

The temperedness of 11^ follows from Proposition 4.3.1 of p^J, or Proposition 14. 91 below. 

Another type of measures in are the stationary empirical fields that are defined 
as follows. Let A C be an open square of side length L, and for lo G let 
^K,pev = {x + Li : x G w, i G Z'^} be its periodic continuation. The associated stationary 
empirical field is then given by 

(2-6) Pa,,^ ^ W\ Ia '^''""'^'P'^' ■ 

It is clear that Pa,i^ is stationary. In addition, it is tempered because 2q{t) < diamA 
for each triangle r G D(a;A^per)- 



4 



2.3 The topology of local convergence 

In contrast to the traditional weak topology on the set of stationary point processes, 
we exploit here a finer topology, which is such that the intensity is a continuous function, 
but nonetheless the entropy density has compact level sets. 

Let £ denote the class of all measurable functions / : Q — > M which are local and 
tame, in that there exists some bounded Borel set A C such that f = f o pr^ and 
I /I ^ b{l + A'^^) for some constant b = b{f) < oo. The topology of local convergence 
on ^0 is then defined as the weak* topology induced by C, i.e., as the smallest topology 
for which the mappings P ^ J f dP with f € C are continuous. By the definition of 
the intensity, it is then clear that the mapping P z{P) is continuous. 

A further basic continuity property is the fact that the centred tile distribution /ip 
depends continuously on P. Let Cq be the class of all hounded measurable functions on 
the space Tq of centred tiles, and =5^ the associated weak* topology on the set A^(7o) 
of all finite measures on Tq- (This is sometimes called the r-topology.) 

Proposition 2.2 Relative to the topologies S'c and JTq introduced above, the mapping 
P ^ fip from ^0 to 7W(7o) is continuous. 

This result will be proved in Section 14.11 It takes advantage of the linear complexity of 
finite Delaunay triangulations, and therefore relies on the planarity of our model. 

2.4 The entropy density 

Given a point process P £ and a bounded Borel set A in M^, we let I\{P,Il^) 
denote the relative entropy (or Kullback-Leibler information) of Pa relative to H^. By 
the independence properties of 11^, these quantities are subadditive in A, which implies 
that the limit 

(2.7) Iz{P)= lim /A(P;n^)/|A| G [0,oo] 

A|^oo 

exists and is equal to the supremum of this sequence. For our purposes, it is sufficient to 
take this limit along a fixed sequence of squares; for example, one can take squares with 
vertex coordinates in Z + 1/2. The claim then follows from the well-known analogous 
result for lattice models [101 Chapter 15] by dividing into unit squares. Iz is called 
the (negative) entropy density with reference measure 11^. 

We set I = Ii. Each 1^ differs from I only by a constant and a multiple of the 
particle density. In fact, an easy computation shows that 

(2.8) I{P) = h{P) + 1 - z + z{P) log z for all z > and P £ S^q . 

A crucial fact we need later is the following result obtained in Lemma 5.2 of |14j . 

Lemma 2.3 In the topology each Iz is lower semicontinuous with compact level 
sets {Iz < c}, c > 0. 
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2.5 Triangle interactions 



This paper is concerned with point processes with a particle interaction which is induced 
by the associated Delaunay triangulation. We stick here to the simplest kind of interac- 
tion, which depends only on the triangles that occur in each configuration. Specifically, 
let (/3 : To — > M be an arbitrary measurable function. It can be extended to a unique 
shift-invariant measurable function (p on T via V'(t) := (/?(t — c(r)), t £ T. Such a (p 
will be called a triangle potential We will assume throughout that (p is bounded, in 
that 



for some constant c<^ < oo. In Theorem 13.41 we will need the following additional 
condition to prove the temperedness of Gibbs measures. Let us say that a triangle 
potential if is eventually increasing if there exist a constant < oo and a measurable 
nondecreasing function ip : [r^, cxd[ M such that <^(r) = iI){q{t)) when q{t) > r^. This 
condition is clearly satisfied when is constant for all triangles r with sufficiently large 
radius ^>(t). 

Example 2.4 Here are some examples of triangle potentials. For each triangle r S T 
let 6(r) = I YIxgt ^ be the barycentre and A{t) the area of r. Examples of bounded 
(and scale invariant) interactions that favour equilateral Delaunay triangles are 



with /3 > 0. Of course, many variants are possible; e.g., one can replace the barycentre by 
the centre of the inscribed circle. By way of contrast, to penalise regular configurations 
one can replace the ipi^s by their negative. 

The triangle potentials ipi above are not eventually increasing. But each triangle 
potential ip can be modified to exhibit this property by setting 



with r > and K a suitable constant. When K is large, one has the additional effect 
of favouring small circumcircles. 

Remark 2.5 The type of interaction introduced above is the simplest possible that is 
adapted to the Delaunay structure. In particular, we avoid here any explicit interaction 
^l' along the Delaunay edges. This has two reasons: First, we might add a term of 
the form | X]eCT-carde=2 fedgei^) to the triangle interaction ip. Such a term would take 
account for an edge interaction ipedge whenever D(a;) is a triangulation of the full plane. 
Secondly, we often need to control the interaction over large distances; the condition of (p 
being eventually increasing is tailored for this purpose. It is then essential to define the 
range in terms of triangles rather than edges. Namely, if a configuration u is augmented 
by a particle at a large distance from oj, the circumcircles of all destroyed triangles must 
be large, but their edges can be arbitrarily short. So, a large-circumcircle assumption 
on the triangle potential allows to control this effect, but a long-edge asumption on an 
edge potential would be useless. 



(2.9) 



99i(t) =/3|c(r)-6(r)|MT) or <^2(t) = -/3 A(t)Mt)2 
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3 Results 



Let (p he a fixed triangle potential. We assume throughout that ip is hounded, see (I2.9p . 
but do not require in general that </? is eventually increasing. For each bounded Borel set 
A C and each configuration w G 17, the associated Hamiltonian in A with boundary 
condition u) is then defined for arbitrary ( £ hy 

(3.1) ^A,..(c)= E '^(^)- 

tGD{CaUwa':): -B(r)nA7^0 

It is always well-defined since the sum is finite. Note that H\^^{^) 7^ in general. 
For defining the associated Gibbs distribution we need to impose a condition on the 
boundary condition lo. 

Definition. Let us say a configuration a; G 17 is admissible if for every bounded Borel 
set A there exists a bounded Borel set h.{oj) D A such that B{t) C A(a;) whenever Q ^Vt 
and T G D(CaUwa=) is such that B{t) n A 7^ 0. We write $7* for the set of all admissible 
configurations. 

In Corollary [12] we will show that = 1 for all P e with P({0}) = 0. Suppose 

now that w G 17* and A is any bounded Borel set. Lemma |2. II then shows that 

where A(a;) is as above. This in turn implies that for each z > the associated partition 
function 

(3.2) Za,,,^ = J e~^^- dUl 

is finite. We can therefore define the Gibbs distribution with activity z > by 

(3.3) Ga,zAA) = ZXi^ I 1a(C U ujac) e-^A-(^) Ui{dC) , AeJ'. 

The measure Ga,z,lj depends measurably on uj and thus defines a probability kernel from 
(17*, J-aO to (17, J^). 

Definition. A probability measure P on (17,.?-") is called a Gibbs point process for the 
Delaunay triangle potential (p and the activity z > 0, or a Delaunay- Gibbs measure for 
short, if P(17*) = 1 and P = / P{du)) Ga,z,uj for all bounded Borel sets A C . We 
write ^0(2;, (^) for the set of all stationary Gibbs measures for and z, and '^^{z,ip) 
for the set of all tempered stationary Gibbs measures; recall Eq. ()2.4p . 

The above definition corresponds to the classical concept of a Gibbs measure, which 
is based on the location of points. We note that an alternative concept of Gibbs measure 
that considers the location of Delaunay triangles has been proposed and used by Zessin 
[2T] and Dereudre [7]. 

Intuitively, the interaction of a configuration in A with its boundary condition oj 
reaches not farther than the set A(li;) above, which guarantees some kind of quasilocality. 
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So one can expect that a limit of suitable Gibbs distributions Ga,z,lj as A | should 
be Gibbsian, and the existence problem reduces to the question of whether such limits 
exist. Our approach here is to take the necessary compactness property from Lemma 
\2.'i\ the compactness of level sets of the entropy density. In fact, we even go one step 
further and show that the stationary Gibbs measures are the minimisers of the free 
energy density. Since such minimisers exist by the compactness of level sets, this solves 
in particular the existence problem. The free energy density is defined as follows; recall 
the definition of the centred tile distribution fip before ()2.3p . 

Definition. The energy density of a stationary point process P G ^@ is defined by 
1>(P) = fy^dfip = \A\'^ j P{duj) ' 

•' •' rGD(a;):c(T)eA 

where A is an arbitrary bounded Borel set. The free energy density of P relative to 
is given by Iz{P) + ^{P)- 

The definition of $ is justified by Proposition 13.61 below. Here are some crucial facts 
on the free energy density, which will be proved in Subsection 14.11 

Proposition 3.1 Relative to the topology on ^' is continuous, and each + $ 
is lower semicontinuous with compact level sets. In particular, the set ^q{z,^p) of all 
minimisers of Iz + ^ is a non-empty convex compact set, and in fact a face of the simplex 

Next we observe that the elements of ^q{z,(p) are nontrivial, in that the empty 
configuration G has zero probability; this result will also be proved in Subsection l4.1i 

Proposition 3.2 For all z > we have ^0 ^ ^q{z,{p), and thus -P({0}) = for all 
P G ^e{z, 

Our main result is the following variational characterisation of Gibbs measures. 

Theorem 3.3 Let ip he a hounded triangle potential and let z > 0. Then every min- 
imiser of the free energy density is a stationary Gihbs measure. That is, the identity 
.J^Q{z,ip) (Z'^Q,{z,ip) holds. In particular, Gibbs measures exist. Conversely, every tem- 
pered stationary Gihhs measure is a minimiser of the free energy density, which means 
that ^q{z, ip) C .Mq^z, Lp). 

The proof will be given in Subsections 14.21 and 14.51 Theorem 13.31 raises the problem 
of whether all stationary Gibbs measures are tempered. It is natural to expect that 
?^e(-2^) V') = ^e'(-^' '^^^ succeed to prove this in general. In fact, we even 
do not know whether 'S'^{z,p>) is always non-empty. But we can offer the following 
sufficient condition, which will be proved in Subsection 14.61 

Theorem 3.4 Suppose ip is eventually increasing and let z > 0. Then every stationary 
Gihhs measure is tempered, so that ^q{z,p) = ^q{z,p). 

Combining Theorems 13.31 and 13.41 we arrive at the following result. 
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Corollary 3.5 Suppose (p is bounded and eventually increasing, and let z > 0. Then 
the minimisers of the free energy density are precisely the stationary tempered Gihhs 
measures. That is, ^q{z, (p) = '^q{z, (p) = '^^q'{z, (p) for all z > 0. 

The proof of Theorem [331 is based on an analysis of the mean energy and the pressure 
in the infinite volume limit when A | M^. For simplicity, we take this limit through a 
fixed reference sequence, namely the sequence 

(3.4) A„ = ] -n-i,n + i[' 

of open centred squares. We shall often write n when we refer to A„. That is, we 
set ujn = t^A„, Pn = Pa„, Rn,uj = -Ra„,i^, Hn,uj = -f^A„,a;, and SO on. We also write 
Vn = |A„| = (2n + 1)^ for the Lebesgue measure of A„. Our first result justifies the 
above definition of $(-P). Besides the Hamiltonian ()3.ip with configurational boundary 
condition iv, we will also consider the Hamiltonian with periodic boundary condition, 
namely 

(3.5) Hn,per{u;) := Vn^Rn,,,) = Yl '^(^)- 

rgD(a;„,pcr): c(T)eA„ 

By definition, we have Hn,-per{^) = 0. Applying Lemma [2TT] and using (j2.9p . we see that 
|-f^n,per| < VnC^p 2 z{Rn) = 2c^ Nn- The following result will be proved in Subsection 14.11 

Proposition 3.6 For every P E we have 

lim [ Hn,pcrdP = '^{P). 

n—foo J 

Moreover, if P is tempered then 

lim [ Hn,MP{duj)=^P). 

n— >oo J 

Finally we turn to the pressure. Let 

Z — p-Hn,'pcr ^TT2 

be the partition function in A„ with periodic boundary condition. 
Proposition 3.7 For each z > 0, the pressure 

p{z, if) := lim log Zn,z,per 

n—>oo 

exists and satisfies 

(3.6) p{z, ^)=- mm [h{P) + $(P)] . 

Proof: This is a direct consequence of Theorem 3.1 of |r4j because ^> is continuous by 
Proposition 13.11 

A counterpart for the partition functions with configurational boundary conditions 
follows later in Proposition 14. Ill Let us conclude with some remarks on extensions and 
further results. 
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Remark 3.8 Large deviations. The following large deviation principle is valid. For 
every measurable A C 

limsupw;;;^ logG„,^,per(i?n ^ A) < -m.iIz^^{c\A) 

n— >oo 

and 

liminf log Gn.z,pcT{Rn ^ A) > - inf /^^^.(int A) , 

n— >oo 

where Gn,z,per = p^j.e~^"-'P''^Il^ is the Gibbs distribution in A„ with periodic bound- 
ary condition, 1^^^ = + $ + p{z, is the excess free energy density, and the closure 
cl and the interior int are taken in the topology =3^. Since ip is bounded so that $ is 
continuous, this is a direct consequence of Theorem 3.1 of [14j . 

Remark 3.9 Marked particles. Our results can be extended to the case of point par- 
ticles with marks, that is, with internal degrees of freedom. Let E be any separable 
metric space, which is equipped with its Borel cr-algebra and a reference measure u, and 
ft the set of all pairs (D = {uj,auj) with a; G and a^j G E'^ ■ In place of the reference 
Poisson point process 11^, one takes the Poisson point process 11 on with intensity 
measure zX u, where A is Lebesgue measure on M?. For ci) G let 

D(w) = {f = (r, ar) -.T G D(w), cr^ = |^} . 

Of course, the centre, radius and circumscribed disc of a marked triangle r are still 
defined in terms of the underlying r. In the germ-grain representation, Tq is replaced 
by the set Tq of all centred r. The tile distribution up of a stationary point process P on 
Q is a finite measure on Tq and is defined by placing bars in (j2.3p . A triangle potential 
is a bounded function if on Tq. Such a (/9 is eventually increasing if ip{f) = iI)[q{t)) 
for some nondecreasing ijj when q{t) is large enough. It is then easily seen that all our 
arguments carry over to this setting without change. 

Remark 3.10 Particles with hard core. There is some interest in the case when the 
particles are required to have at least some distance ro > 0. This is expressed by adding 
to the Hamiltonian (jS.ip a hard-core pair interaction term -^^^^^(C) which is equal to 
oo if |x — y| < ro for a pair {x, y} C Ca U uj\c with {x, y} n A / 0, and zero otherwise. 
Equivalently, one can replace the configuration space by the space 

X — y\ > ro for any two distinct j;, y G w} 

of all hard-core configurations. The free energy functional on then takes the form 
F^" := + $ + ^>^^ where 

^''"'{P) = oo : < |x| < ro for some x G w) = oo \ 

for P G with Palm measure P^; here we use the convention ooO = 0. We claim 
that our results can also be adapted to this setting. In particular, the minimisers of F^'^ 
are Gibbsian for z and the combined triangle and hard-core pair interaction, and the 
tempered Gibbs measures for this interaction minimise F^'^. We will comment on the 
necessary modifications in Remarks 14.41 and 14.121 



10 



Combining the extensions in the last two remarks we can include the following 
example of phase transition. 



Example 3.11 The Delaunay-Potts hard-core model for particles with q > 2 colours. 
In the setup of Remark 13.91 we have E = {1, . . . ,q}, and the triangle potential is 



where /? > is the inverse temperature and ri > is an arbitrary interaction radius. 
If one adds a hard-core pair interaction with range tq < ri/^/2 as in Remark 13. 101 this 
model is similar to the model considered in (Instead of a triangle potential, these 
authors consider an edge potential along the Delaunay edges that do not belong to a tile 
T of radius ^(r) > ri.) Using a random cluster representation of the triangle interaction 
as in [15] and replacing edge percolation by tile percolation one finds that the methods 
of [1] can be adapted to the present model. Consequently, if z and f3 are sufficiently 
large, then the simplex ^q^z, (p) = l^e(-2, V') ^las at least q distinct extreme points. 

4 Proofs 

4.1 Energy and free energy 

We begin with the proof of Proposition[221 which states that the centred tile distribution 
/ip depends continuously on P. The continuity of the energy density $ and the lower 
semicontinuity of the free energy density + <& then follow immediately. 

Proof of Proposition 1 2. Let {Pa) be a net in that converges to some P G 

We need to show that J g d^p^ ^ J g dfip for all g £ Cq. We can assume without loss 

of generality that < g < 1. 

We first consider the case that g has bounded support, in that g < l{g < r} for 
some r > 0. Let A C be any bounded set of Lebesgue measure |A| = 1. Also, let 



In view of (j2.3p , we then have / g d^Q = J f dQ for all Q G , and in particular for 
Q = Pa and Q = P. By the bounded support property of g, f depends only on the 
configuration in the r-neigbourhood A^' := {y £ M."^ : \y — x\ < r for some x € A} of 
A. That is, / is measurable with respect to J^/\r. Moreover, / < 2N/\r by Lemma l2.H 
so that f £ C. In the present case, the result thus follows from the definition of the 
topology ^c- 

If g fails to be of bounded support, we can proceed as follows. Let e > be given 
and r > be so large that fJ.p{Q > r) < e. Since ||/Up^|| = 2z{Pa) — > 2z{P) = 
||/ip|| and fJ-p^ig < t) ^ fJ'p{6 ^ ''') by the argument above, we have fJ-p^ig > r) — > 
fip{g > r). We can therefore assume without loss of generality that ^p^ig > r) < e for 




(3 if g{T) < ri and a-r is not constant, 
otherwise. 



/M= Yl 5(r-c(r))lA(c(r)). 
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all a. Using again the first part of this proof, we can thus write 
y g dfip - e < j g l{g<r} d^ip = lim j g l{g<r} dfip^ 
< liminf / gd^p^ < limsup / gdfip^ 

° y a J 

- ^™ / ^ l{e<r} d^ip^ + e = j g l{g<r} d^ip + e < j g dfip + e . 

Since e was chosen arbitrarily, the result follows. 

We now turn to the properties of the free energy density. 

Proof of Proposition ] 3. 1{ As ip belongs to Cq, the continuity of ^ follows immediately 
from Proposition 12.21 By Lemma [2.3l we can also conclude that Iz + ^ is lower semicon- 
tiuous. Moreover, hypothesis (j2.9p implies that the level set {/^ + <I> < c} is contained 
in {Iz < c + 2cipz{-)}, which by p. 81) coincides with the compact set {I^/ < c + z' — 1} 
for z' = z exp(2c^). 

Let P = (50 G ^0 be the Dirac measure at the empty configuration. Then fip = 
and thus = 0. On the other hand, Iz{P) = z. This means that + <^ is not 

identically equal to +oo on S^q and thus, by the compactness of its level sets, attains 
its infimum. To see that the minimisers form a face of it is sufficient to note that 
/j; + <l> is measure affine; cf. Theorem (15.20) of [lOj. 

Next we show that the minimisers of the free energy are nondegenerate. 

Proof of Proposition [37B: The second statement follows from the first because Mq{z, ip) 
is a face of ^e- For, suppose there exists some P £ M.Q{z,ip) with P({0}) > 0. 
Then 60 appears in the ergodic decomposition of P, which would only be possible if 
(50 G Me{z,(p). 

To prove the first statement we note that $(11") < C(p||^n"|| = 2c^u for all n > 0. 
Therefore, if z > is given and u is small enough then 

(4.1) /,(n") + $(n^) <z-u + ulog{u/z) + 2c^u <z = Iz{50) + $((50) , 
so that (5^ is no minimiser of the free energy. <0 

Finally we show that the energy density $ is the infinite volume limit of the mean 
energy per volume. 

Proof of Proposition \3.(k We begin with the case of periodic boundary conditions. 
For every P e =^e, we have j Hn,pcrdP = J ^{Rn)dP = $(Pi?„). It is easy to 
see that PRn — > P, cf. Remark 2.4 of [H]. Since $ is continuous, it follows that 

^PRn) ^ HP). 



12 



Next we consider the case of configurational boundary conditions and suppose that 
P is tempered. Applying p.3p we obtain for each n 

/ P{dw)Hn,u;{uj) = I P{dw) ^ (^(r - c(r))l{B(r-c(r))n(A„-c(r))5^0} 
(4.2) = j Hp{dT) (p{t) j dx l{ij(r)n(A„-x)5^0} 

where A.^^^^ is the ^?(T)-neigbourhood of A„. Now, for each r we have 



|A^W|K = l+Mr)/V^+TTQ{Tf/Vn ^ 1 



as n — > oo. 



In view of (12. 4p and ()2.9p , we can apply the dominated convergence theorem to conclude 
that 



^P) = lim [ Hn,U^)P{du;), 

n— »oo J 



as desired. 



4.2 The variational principle: first part 

In this section we shall prove that each minimiser of the free energy is a Delaunay-Gibbs 
measure. We start with an auxiliary result on the 'range of influence' of the boundary 
condition on the events within a bounded set. Let A C be a bounded Borel set and 
a; G ri. Writing B{x,r) for the open disc in with center x and radius r, we define 

7^A(a;) = {r > : card WB{x,r)\A > 1 for all x G s.t. B{x, r) n A 7^ 0} 

and rj\{uj) = miTZ/\{u;), where inf0 := 00. Let A^*" = [J^^^ B{x,2r) be the open 
2r-neigbourhood of A. We then observe the following. 

Lemma 4.1 (a) For all r > 0, {ta < r} G J^^2r\/^. In particular, fa is Ti^<^- 
measurable. 

(b) For all P e we have P{rA = 00) = P({0}). 

Proof: (a) Let 7^a('^) be defined as 7^a('^)) except that the discs are required to have 
rational centres x £ Q^. Then TZ/\{uj) C 72.a(w). Moreover, if r < r' then every open 
r'-disc intersecting A contains an r-disc with rational center and intersecting A, so that 
r' £ TZ/\{uj) when r G 7?.a(w). This shows that 

{rA<r}= IJ fl {A^ij(.,s)\A > 1} , 

seQ: 0<s<r xeQ^: B(x,s)nA^I!) 

and the last set certainly belongs to .FA2r\^A- 

(b) Since {0} = flreNi-^sCo.r) = 0} C {fa = cxd}, it is sufficient to prove that 
{ta = CXD, A''^(o,r) ^ 1} has measure zero for all r > and P G ^e- However, if 
{ta = cxd} occurs then there exists a cone C with apex at some point in the closure of 
A and an axis in one of finitely many prescribed directions such that Nq = 0, while 
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the Poincare recurrence theorem imphes that Nq = oo for each such C almost surely 
on {NB(o^r) ^ !}• This contradiction gives the desired result. 

As an immediate consequence we obtain that each nondegenerate stationary point 
process is concentrated on the set ri* of admissible configurations. 

Corollary 4.2 The set fi* of admissible configurations is measurable (in fact, shift 
invariant and tail measurable), and P(0,*) = 1 for all P G with P({0}) = 0. 

Proof: This is immediate from Lemma WA\ because VL* = nn>i{'^A„ < oo}. 

Next we state a consequence of Proposition 13.71 

Corollary 4.3 For every P G ^Q{z,ip), we have 

lim V~^In{P]Gn,z,per) = . 

n— >oo 

Proof: By the definition of relative entropy, 

In{P\Gn,z,pcT) = In{P]^'^) + j Hn,pcT dP + log Zn,z,peT ■ 

Together with Propositions 13.6 1 and 13 .71 this gives the result. 

We are now ready to show that the minimisers of + $ are Gibbsian. 

Proof of Theorem \3.3\, first part: We follow the well-known scheme of Preston [19] (in 
the variant used in p3], Section 7). Let P £ ^q{z,(p), / be a bounded local function, 
A a bounded Borel set, and 



/aH = j f{C)GA,z,UdO, <^^^- 



We need to show that J f dP = / /a dP. Let rA be the range function defined above, 
and for each r > let lA,r = l{rA < ''^j and A^^ be the 2r-neigbourhood of A. By 
Lemma im (a). lA,r is measurable with respect to .F^2r\^. Moreover, if rA(a;) < r then 
Ha,u){C,) = HA,Lu^2r\^{0 for all Q G f^A- So, if r is so large that / is J^^^^'-measurable, 
we can conclude that lA,r /a is .^Aar^^-measurable. 

Now we apply Corollarv 14.31 which states that lim^^oo v~^Ia„{P', Gn,z,per) = 0. By 
shift invariance, this implies that Pa G'A,2,per with a density g\ for each sufficiently 
large square A. In particular, for any A' C A we have Pa' ^ (G^A,z,per)A' wih density 
gA,A'i^) ~ I G A\A',z,ujn A' (dC) 9 a{0- Corollary 14.31 implies further that for each 6 > 
there exists a square A and a Borel set A' with A^*" C A' C A such that f \gA,A' — 
5a,A'\aI dG\^z,per < S', cf. Lemma 7.5 of [I3]. Now we consider the difference 

j lA,r (/ - /a) '^^ ^ j ^\r{9A,A' f ~ 9A,A'\A /a) dGA,z,per ■ 

Since GA,z,pcr = / GA,z,pcr{du!) Ga,z,ui and lA.rffA.A'XA is .?^A\A-™easurable, we can 
conclude that 



J lA,r 5A,A'\A /a '^G'A,z,pcr = y ^A,r 9 A,A'\A f dG A, 



2, per • 
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By the choice of A and A', we can replace the density gA,A'\A the last expression by 
5a,A' making an error of at most 6. We thus find that J lA,r (/ — /a) dP < 5. Letting 
(5 — > and r — > oo, we finally obtain by the dominated convergence theorem that 

/ if-fA)dP = 0. 

J {rA<oo} 

This completes the proof because P(rA = oo) = P({0}) = by Lemma 14.1( b) and 
Proposition 13. 2i 

Remark 4.4 Here are some comments on the necessary modifications in the hard-core 
setup of Remark 13.101 In analogy to Proposition 13. 7^ one needs that 

lim y-^ log / e~^"'P<='~-^""p- dffi = - min \lJP) + ^(P) + '^'^%P)] . 

n^oo J PG.^e '- 

This follows directly from Propositions 4.1 and 5.4 of [12] because $ is continuous. 
Corollary 14.31 therefore still holds for the periodic Gibbs distributions with additional 
hard-core pair interaction. One also needs to modify the proof of Proposition 13.21 in 
that the Poisson processes 11" should be replaced by the Gibbs measure with activity 
u and pure hard-core interaction. P" is defined as the limit of the Gibbs distributions 
^n%,per ^'^^ the periodic hard-core Hamiltonians -ff^^per- By Proposition 7.4 of [12], P" 
exists and satisfies 

/„(P") = - hm v-^\og [ e-^"^p- dUl < - lim logn;^({0}) =u. 

n—>oo J n—>QO 

Together with ^M) we find that 

4(P") + cD(P") < z + z(P") [log(u/z) + 2cJ , 

which is strictly less than z when u is small enough. Since <I)'^'^(P") = ^^^[Sd)) = 0, it 
follows that the minimisers of Iz + ^ + ^^'^ are non-degenerate. No further changes are 
required for the proof of the first part of the variational principle. 



4.3 Boundary estimates 

We now work towards a proof of the reverse part of the variational principle. In this 
section, we control the boundary effects that determine the difference of Hn^per and 
Hn^uj- The resulting estimates will be crucial for the proof of Proposition 14.111 For 
every io G Q and every Borel set A let 

Sa{uj) = {r G D{uj) : P(r) n A / and B{t) \ A / 0} 

be the set of all triangles r E D(w) for which P(t) crosses the boundary of A. We start 
with a lemma that controls the influence on Sa when two configurations are pasted 
together. 

Lemma 4.5 Let A be a (not necessarily bounded) Borel set in M?, C G 17* U {0} a 
configuraton with (qa = 0, and w G 17. Then for each r G Sa(Ca U ua") and each 
X G TA there exists some t' G Sa(C) '^^^^ x G r'. 
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Proof: Let A, and uj be given. If ( is empty, there exists no x G C Caj so that the 
statement is trivially true. So let ( G Q* and suppose there exists some r E Sa(CaUwa<:) 
with TA / 0. Let x G ta- Since (q^ = 0, x does in fact belong to the interior of A. 
This implies that B{t') n A 7^ for each r' € D((^) containing x. Therefore we only 
need to show that B{t') \ A 7^ for at least one such r'. Suppose the contrary. Then 
B{t') C a whenever x € r' € D(C). This means that the Delaunay triangles containing 
X are completely determined by Ca- This gives the contradiction 

/ {r G Sa(Ca U waO : r 9 x} = {r' G Sa(C) : r' 9 x} = , 

and the proof is complete. 

The following proposition is the fundamental boundary estimate. It bounds the 
difference of Hamiltonians with periodic and configurational boundary conditions in 
terms of := cardSA^. 

Proposition 4.6 There exists a universal constant 7 < 00 such that 

\Hn,per{0 " Hn,UC)\ < 7C^ (^n(^) + 5„(C)) 

for all n>l and all C,cj G il* U {0} with C,d\„ = ^d\n — 0- 

Proof: Let n, ^, u; be fixed and 

A = {r G D(Ca„ U l^ac ) : B{t) n A„ / 0} , B = {r G D(C„,pcr) : c(r) G A„} . 

In view of (f3T]) and ([33]) we have Hn^uiC) = Y^reA'^i'^) ^^'^ -f^n,pcr(C) = Z^reB'^('^)- 
Since if is bounded by c^, we only need to estimate the cardinalities of A \ B and B \ A. 
We note that A\B C SA„(CA„Ua;A= ) and B\A C SA„(Cn,per)- So we can apply Lemma [13] 
to both A = An and A = to obtain that the set of points belonging to a triangle in 
A \ B is contained in the set of points belonging to a triangle of Sa„(C) U Sa„(w). Hence, 
card ( IJ^g^^^g r) < 3(S'„(C) + 5„(c<j)). By Lemma \2.1\ it follows that card (A \ B) < 

6(5„H + s„(C)). 

To estimate the cardinality of B \ A we may assume that Cn 7^ 0- The periodic con- 
tinuation Cn,per then contains a lattice, and this implies that every triangle of D(^„^per) 
has a circumscribed disc of diameter at most \/2vn- Hence, each r G B \ A is con- 
tained in A5„+2, the union of 5^ translates of A„ (up to their boundaries). Applying 
Lemma 14.51 to each of these translates we conclude that the number of points that be- 
long to a triangle of B \ A is bounded by 3 • 5^ 5'n(C)- Using Lemma |2. II again we find 
that card(i3 \ ^) ^ 150(S'„(a;) -|- Sn{C))i ai^cl the result follows with 7 = 156. 

The following immediate corollary will be needed in the proof of Theorem 13.31 

Corollary 4.7 There exists a constant C < 00 such that \Hn,uj{9)\ < C Sn{(^) for all 
n > 1 and uj G Q* with wqa^ = 0- 

The next proposition exhibits the fundamental role of the temperedness condi- 
tion (|2.4p combined with stationarity for controlling the boundary effects. 

Proposition 4.8 For every P G r^pQ, v~^Sn in L^{P) and P-almost surely. 
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Proof: For each z G we consider the shifted unit square C{i) = Aq + i and define the 
random variable 

Zi = cardjr £ D(-) : B{t) f] C{i) / 0}. 

Then 

(4.3) Sn< Yl 

i£ln\In-l 

where /„ = A„ n (Z^ + (^) |))- Note that card In = Vn- As in (|4.2p we have 

J ZodP = j fIp{dT) |Ag(")| = j f,p{dT) {l+MT)+7T0{Tf). 

The last term is finite by the temperedness of P. So, each Zi is P-integrable. Since 
Zi = Z^odi, the two-dimensional ergodic theorem implies that SiG/„ converges to 
a finite limit Z, both P-almost surely and in L^{P). This implies that 'Yli^i^Xin i 
tends to zero P-almost surely and in L^{P). The result thus follows from ()4.3|) . 

4.4 Temperedness and block average approximation 

Our first result in this subsection is a sufficient condition for temperedness in terms of 
vacuum probabilities. For P £ let 

(4.4) Vk{P) = ess sup P{Na, = 0\Tai) 

be the essential supremum of the conditional probability that A^ contains no particle 
given the configuration outside. 

Proposition 4.9 Every P £ satisfying 

(4.5) Y.'^kVkiP) z{P) < ^ 



k>0 

is tempered. 



Proof: We can assume that P ^ Sij^ because otherwise the result is trivial. For each 
/c > 1 we consider the shifted squares A^,(i) = A^ + {2k + i G Z^, as well as the event 

Ak = {iVAfe(j) / for ah i G 1? with ||z||oo < l}. 

Since P 7^ (50, it is clear that P{A\^ — > 1 as A; — > 00. Thus we can write 

y card(r G D(a;) : G B{t)) P{du 



< ^ / card(T G D(w) : G B{T))P{du)) 



with the convention Aq = 0. Now, if w G then each circumscribed disc containing 
of a triangle r G 0{uj) has a diameter not larger than 2y^2t^, so that each such r in 
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fact belongs to D{uj\^i^_^^^). Lemma [2. II thus shows that the number of such r is at most 
2-^7fc+3(<^)- The last sum is therefore not larger than 

k>l "^k-l k>li£.I?-.\\i\\^<l 

In view of the stationarity of P, the last integral is bounded by Vk-i{P) vjk+3 z{P) = 
7^ Vk Vk-i{P) z{P)- So we arrive at the estimate 

/ card(r G D(a;) : G B{t)) P{duj) < 2 • 7^ 3^ ^Ufc Vfc_i(P) z{P) . 
•' k>i 

Together with <\2.b\) and assumption ()4.5p . this implies the temperedness of P because 
ffc ~ as /c ^ oo. <0 

The second result concerns the approximation of stationary measures in terms of 
tempered ergodic measures. This approximation uses the block average construction 
first introduced by Parthasarathy for proving that the ergodic measures are dense in 
^e; cf. [IQi Theorem (14.12)], for example. 

Proposition 4.10 Let z > and Q e be such that Iz{Q) + ^{Q) < oo. Then for 
each e > there exists some tempered Q-ergodic Q £ such that Iz{Q) < Iz{Q) + ^ 
and $(Q) < ^{Q)+e. 

Proof: Let Q E =^^e be given. We can assume that Q ^% because otherwise we can 
choose Q = Q. For n > 1 let Q"*^ denote the probability measure on Q relative to which 
the particle configurations in the blocks A„ + (2n + i G Z^, are independent with 
identical distribution Qn = Q o pr^^. (In particular, this means that the boundaries of 
these blocks contain no particles.) Consider the spatial average 



J An 



It is clear that Qji*^"""^ G ^e- It is also well-known that Qjj'^"''^ is 0-ergodic; cf. [IHl 
Theorem (14.12)], for example. By an analogue of [101 Proposition (16.34)] or [14^ 
Lemma 5.5], its entropy density satisfies 

(4.6) h{Q-^-n = Vn'HQ;K)- 

So, IziQn'^'^'') < Iz{Q) + e when n is large enough. Moreover, the same argument as in 
[m Lemma 5.7] shows that Q^^~^^ ^ Q in the topology ■S/'c- By Proposition 12.21 <I> is 
continuous. We thus conclude that ^{Q'^f^'^^) ^{Q), whence ^{Q"^^'''"') < $(P) +e 
for large n. It remains to prove that each Q^^'^"^ is tempered. Let n be fixed and 
k > £(2n+l) for some i > 1. We claim that Vk{Q'^'^~''^) < q'"'-' with q = Q{Nn = 0). 
Indeed, for each x G A„ we have A^, + x D A„_,_(£_;^)(2n+i)i ^iid the latter set consists of 
Vi-i = {21 — 1)^ distinct blocks as above. Letting g be any nonnegative .Fa= -measurable 
function and using the independence of block configurations, we thus conclude that 



A„ 



iid-av 
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which proves the claim. Now, we have q < 1 because Q ^ Sijf. It follows that 

k>2n e>l e{2n+l)<k<{e+l){2n+l) £>0 

for some constant C„ < cxd. Together with Proposition 14.91 this gives the temperedness 
of Qj^'l-'^^. 

4.5 The variational principle: second part 

In this section we will complete the proof of the variational principle. The essential in- 
gredient is the following counterpart of Proposition 13. 71 involving configurational instead 
of periodic boundary conditions. We only state the lower bound we need. 

Proposition 4.11 For every P £ and P -almost every uj we have 

liminf log Z„ ^ ^ > p{z, (/?) . 

71— ^OO 

Proof: By (|3.6p and Lemma I4.1UI it is sufficient to show that 

liminfT;-ilogZ„,,,^ > -h{Q) - ^{Q) 

n—>oo 

for every ergodic Q G ^q". We can assume without loss that the right-hand side is 
finite. Now, since Iz{Q) is finite, Q is locally absolutely continuous with repect to li^ . 
We fix some e > 0, let /„ = dQn/dH^, and consider for every lj £ ^l* the set 

An,cu = {\Hn,cu - Hn,per\/Vn < £, '^{Rn) < $(Q) + £, V'Hog fn < hiQ) + e} . 

Then for sufficiently large n we have 

> g-fn [/4Q)+-I>{Q)+3e] Q(^^ ^) _ 

It is therefore sufficient to show that for P-almost u), (5(^n,a;) ^ 1 as n ^ cxd. By the 
ergodic theorem, $(-Rn) converges to ^{Q) in Q-probability; cf. Remark 2.4 of |14j . By 
McMillan's theorem [Hi [18], Q{v~^ log fn < Iz{Q) + s) —>■ 1 when n tends to infinity. 
Moreover, Propositions 14.61 and 14.81 imply that, for P-almost all uj, \Hn^i^ — Hn^per\/vn 
converges to in L^{Q). This gives the result. 

We can now show that every tempered stationary Gibbs measure minimises the free 
energy density. 

Proof of Theorem 13.31 second part: We follow the argument of [13], Proposition 7.7. 
Let P G ^0^(93, z). On each J^k„, P is absolutely continuous w.r. to 11^ with density 

5n(C) = / P(d^)^^(C) = / P(d^)e-^-(^)/Z„,.,.. 
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Using Jensen's inequality and the Gibbs property of P we thus find that 



= - j P{duj)Hn,M- j P{duj) logZ„,,,^. 

Next we divide by u„ and let n oo. We know from Proposition 13.61 that 
v~'^ J P{du}) Hn,uj{uj) ^{P)- On the other hand, Corollary [12] implies that 
v^^ logZn^z,ui > —z—Cv~^Sn{uj). Using Propositions 14. 8l and l4.11l together with Fatou's 
Lemma, we thus find that 

lim inf u"^ / P{duj) log Zn,z,uj > p{z,f)- 

Therefore Iz{P) < —^{P) + min[/2 + as required. 

Remark 4.12 In the hard-core setting of Remark 13.101 a slight refinement of Propo- 
sition HTTO] is needed. Namely, under the additional assumption that ^^'^{Q) = one 
needs to achieve that also ^^'^{Q) = 0. To this end we fix an integer k > ro/2 and 
define Q^^'^ in such a way that the particle configurations in the blocks A„ + (2n + 
z G Z^, are independent with identical distribution Qn-k, rather than This means 
that the blocks are separated by corridors of width 2k > ro that contain no particles. 
It follows that = 0, and it is still true that limsup„_o^ /^(Qji'^"^^) < h{Q); 

cf. |12t Lemma 5.1]. We thus obtain the refined Proposition 14.101 as before. 

A similar refinement is required in the proof of Proposition 14.111 One can assume 
that (^^''{P) = and = 0, and in the definition of An. ^ one should introduce an 

empty corridor at the inner boundary of A„ to ensure that H^'^^ = H^'^p^j. = on An^uj, 
see [I2l Proposition 5.4] for details. In the proof of the second part of Theorem 13.31 one 
then only needs to note that ^^'^{P) = when P is a Gibbs measure P for the combined 
triangle and hard-core pair interaction. The proof of Theorem 13.41 carries over to the 
hard-core case without any changes. 

4.6 Temperedness of Gibbs measures 

Here we establish Theorem 13.41 By Proposition 14.91 it is sufficient to show the following. 



Proposition 4.13 Let ip he hounded and eventually increasing, z > 0, and P he any 
stationary Gihhs measure for ip and z. Then there exists a constant C > such that 

(4.7) P{Nk = 0\J'Ai)<Cv,' 

for all k > 1. 

To prove this we need an auxiliary result which states that the radii of all cir- 
cumcircles in the Delaunay tessellation must decrease when a point is added to the 
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Figure 1: D{uj) (solid lines) and D(a;U{x}) \ D(a;) (dashed lines). The difference region 
Ax{uj) is shaded in grey. 

configuration. Specifically, let w G il* and x G \ a; be such that lo U {x} is in general 
circular position and x is not collinear with two points of u. We consider the sets 

Q(u;) := D(a;) \ D(w U {x}) = {r G D(w) : B{t) B x} 

and 

C+(w) := D{uj U {x}) \ D{uj) = {tG D{uj U {x}) :tBx}. 
If (r) denotes the convex hull of a triangle r, 

(4.8) A.(^):= U ^^)= U 

is the region on which the triangulations D{uj) and D{uj U {x}) differ; see Fig. [TJ Up to 
the point x, the interior A°(w) of Aa;(a;) is covered by the discs B{t) with r G C+(a;), 
which by definition are free of particles. Consequently, A°(a;) contains no particle of lo, 
so that the vertices of each r G Cx{uj) belong to the boundary dAx{co)- Next, Lemma 
12.11 shows that 

(4.9) cardC+(a;) = cardQ.(Lj) + 2 , 
and for every A B x we have 

(4.10) Ha,^{co U {x}) - Ha^U^) = ^ ip{T)- ^ (^(r). 

TGC+(a;) reCa^iui) 

Here is the monotonicity result announced above. 

Lemma 4.14 Under the conditions above, there exist a subset C'^{uj) C Cx{uj) with 
card(C2^(w) \ C^(w)) < 4 and an injection I from Cx{lo) to C^iuj) such that 

(4.11) q{I{t))<0{t) for allT eC'xiuj). 
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We postpone the proof of this lemma until the end, coming first to its use. 



Proof of Proposition \4.13 : By assumption, ip is eventually increasing. So there exists 



some r^ < oo and a nondecreasing function ■0 such that (^(r) = il;{g{T)) when q(t) > r^. 
Combining Lemma 14.141 and Equations ()4.9p and ()4.10p we thus find that 

(4.12) Hk,^{LO U {x}) < HkA^) + 10 

for all Lij S ri*. A; > 1, and Lebesgue-almost all x G \t(; that have at least the distance 
2r^ from all points of u. Next, let P G ^e(z, (/?). By definition, 

P{Nk = 0|.FA=J(a;) = Z-i^^e— e-^'^-W 

for all w G il*. Let A^.^^ = {(x,y) G A^_2^^ : Ix-y] > 2r^}. Applying (|4.12p twice (viz. 
to and X as well as uj\c u {x} and y) and recalling ()3.2|) we find that 



> |A^^^| e"-""*^ e"-'^'='^(^)"^°'^'^/2 . 



Since |A^^^| ~ as fc — > cxd, the result follows. 

Finally we turn to the proof of Lemma I4.14[ 

Proof of Lemma \4-14\ ' Let Tx be the unique triangle of Cx{uj) containing x in its 
interior, and C^^(u;) the set of all r G C^iuj) that have an acute or right angle at x. 
Note that card(Cj^(tij) \ C^^(a;)) < 3 because the angles at x of all r G C|^(a;) add up 
to 360 degrees. We will associate to each triangle r G Cx{uj) a triangle /(r) G C^{lj), 
except possibly when r = r^; or the candidate for /(r) does not belong to Cx^{u)). Our 
definition of /(r) depends on the number k = kij) of edges e C r with (e) C dIS.x{uj). 
Let Cif^(a') be the set of all r G Cx{'^) that have fc such edges. Since ci^'*(a;) = except 
when Cx{'^) = {''"xji we only need to consider the three cases fc = 0, 1,2. 

The cases fc = 1 and 2 are easy: For every r G C^x\oj) there exists a unique edge e(T) 
such that e(T) U {x} G C+(a;). If in fact e(r) U {x} G C+^(w) we set /(r) = e(r) U {x}; 
otherwise we leave I(r) undefined. Likewise, every r G Ci^ (w) has two edges ei(T) and 
62 (r) in dl^x{^) (in clockwise order, say) and can be mapped to I(r) = ei(T) U {x}, 
provided this triangle belongs to C|^'^(w). The resulting mapping / is clearly injective. 
Moreover, r and /(r) have the edge e(T) (resp. ei(T)) in common, and x G -B(t) 
because r G Ca;(a'). Since /(r) G C+^(w) whenever it is defined, we can conclude that 

e{i{T)) < q{t). 

The case fc = is more complicated because the tiles r G ci'''*(a;) are not naturally 
associated to a tile of Cj^(a;). To circumvent this difficulty we define an injection / from 
C^x\'^) \ {tx} to C^x\uj) such that q{I{t)) < q{t). Each triangle r G C^x\oj) different 
from Tx can then be mapped to the triangle /(r) = e2{I{T)) U {x}, provided the latter 
belongs to C^'^(a;); otherwise /(r) remains undefined. This completes the construction 
of /. (Note that Tx does not necessarily belong to C^x\i^)- However, if it does we have 
no useful definition of I{tx)-) 
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Figure 2: The set C^(u;) for the configuration u of Fig. 1, with a tile r* S Ci: (cj) (hght 
grey), its circumcircle B(t^) (dashed), the associated edges and regions Wi (dark 
grey), and two triangles Tj G C2;(u;) with Tj C Wi with their circumcircles (solid). The 
construction in the proof gives /(r*) = T2. 

To construct / we turn Cx{u!) into the vertex set of a graph Gxito) by saying that 

(2) 

two tiles are adjacent if they share an edge. The set (w) then coincides with the 
set of all leaves of Gx{to), and C^x\uj) is the set of all triple points of Gx{uj). Consider 
a fixed G C^x\'-^) \ {tx}- Since r^< C dAx{uj), the set Ax{uj) \ (r*) splits into three 
connected components. Let Wi = Wj(r^<,x,u;) be the closure of the ith component, 
i = 1,2,3. Any two of these sets intersect at a point of r*, and one of them contains 
X because t^: / r^. Suppose x E W3. For i = 1, 2 let Cj = n be the edge of t^: 
that separates Wi from the rest of Ax{uj); see Fig. [2l We claim that there exists some 
i = i{T^) G {1, 2} such that 

(4.13) £>(t) < g{n) for all r € Cx{oj) with t C Wi . 

The image /(n) of r* can then be defined as the leaf of G^(cl') in with the largest 

circumcircle. (The largest circumcircle condition takes account of the fact that the path 
in Gx{i^) from to /(r*) might contain further triple points.) 

It remains to prove (|4.13p . Since / t^, there exists at least one i such that the 
triangle {x} U Cj has an acute angle at x. We fix such an i and consider any r € Cx{co) 
with T C Wi. There exists at least one point zq £ t that is not contained in the closed 
disc B{t^). Since x S B{t) and (r) is covered by the tiles (r') for r' € C^(ll') with 
(r') n (cj) / 0, we conclude that the line segment s from zq to x is contained in B{t) 
and hits both the circle dB{T^) and the edge (ej). In particular, -B(r) n (ej) / 0. Since 
-B(t) contains no points of a;, we deduce further that the circle 9i?(r) hits the edge (cj) 
in precisely two points zi and Z2. By the choice of the angle of the triangle {zi,x, Z2} 
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at X is acute. Since x € B{t), it follows that the angle of the triangle {zi, zq, Z2} at zq is 
obtuse. Consequently, if we consider running points such that yo runs from zq to the 
point sn5i3(r*) and the edge {yi,y2} from {zi,Z2} to Cj, the associated circumcircles 
2/2}) run from S(t) to B{t^), and their radii Q{{yi,yo,y2}) must grow. This 
proves that g{T) < £)(t*), and the proof of (j4.13p and the lemma is complete. {> 

Acknowledgment. We are grateful to Remy Drouilhet who brought us together 
and drew the interest of H.-O. G. to the subject. 
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